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In a previous paper [1], we have discussed the non-perturbative tuning of the chirally 
£SJ rotated Schrodinger functional (xSF). This tuning is required to eliminate bulk 0(a) 

cutoff effects in physical correlation functions. Using our tuning results obtained 
in [1] we perform scaling and universality tests analyzing the residual 0(a) cutoff 
\Q effects of several step-scaling functions and we compute renormalization factors at 

the matching scale. As an example of possible application of the %SF we compute 
the renormalized strange quark mass using large volume data obtained from Wilson 
twisted mass fermions at maximal twist. 
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1 Introduction 



The Schrodinger functional (SF) scheme [2"j3jT4] has been widely employed for performing 
non-perturbative renormalization and scaling studies. For an incomplete list of references 
see 
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It is not straightforward to implement a Schrodinger functional (SF) scheme that retains 
the property of automatic 0(a) improvement [H] for massless Wilson fermions. A solu- 
tion to this problem has been proposed recently by Sint [15J, called the chirally rotated 
Schrodinger functional (%SF). In refs. [T6|17Jll8j preliminary studies of the %SF were pre- 
sented and, recently in a companion paper [1] , we have performed a detailed study of the 
non-perturbative tuning of the %SF in a quenched setup. To retain the property of bulk 
automatic 0(a) improvement, two parameters of the lattice action, one denoted by Zf and 
the other, the usual hopping parameter k, have to be tuned non-perturbatively to their 
critical values. 

In ref. [1] we have defined a tuning procedure and we have found that the tuning of 
Zf and k is numerically feasible with standard techniques. We have found that several 
different tuning conditions lead to fully compatible results for all physical correlation 
functions and we have also numerically checked that in the continuum limit the correct 
boundary conditions, i.e. the boundary conditions needed to have bulk automatic 0(a) 
improvement, are recovered. 

The xSF, being compatible with automatic 0(a) improvement, can hence be used for 
renormalizing bare operators computed with Wilson twisted mass fermions at maxi- 
mal twist [IH] and the computation of operator specific improvement coefficients can 
be avoided. In this paper we want to demonstrate this explicitly for the particular case 
of a twist-2 operator. Additionally, we discuss the application of the %SF scheme for the 
renormalization of quark masses. 

The computation of such quantities allows us to perform a test of the universality of 
the continuum limit of lattice QCD and, as a result, to confirm the correctness of the 
continuum limit of the novel xSF scheme itself. The latter can be achieved by comparing 
the continuum limit values of the computed quantities in the %SF formulation to those 
values obtained using the standard version of the SF. Since, as discussed in ref. [T5fT] . 
the SF and the %SF are the same formulations in the continuum, the final results in the 
continuum limit, as obtained from the two formulations, should be exactly the same given 
a common choice of a renormalization prescription. 

The paper is organized as follows. We discuss the continuum limit of the step-scaling 
function (SSF) of the pseudoscalar density, up, in sect. [2] and of the non-singlet twist-2 
operators, o"o 12 and ao 44 , in sect. [3j The continuum limit values are then compared to 
the ones obtained using the SF, with which, as expected, we find agreement. Next, in 
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sect. H] we show the results of the determination of the renormalization factors of the 
pseudoscalar density, Zp, and the twist-2 operators, Z Ql2 and Z Qii , at hadronic matching 
scales and non-zero lattice spacing. For the latter, results are compared against values 
obtained using the standard formulation of the SF with two different regularizations, non- 
improved Wilson and clover improved Wilson fermions. Eventually in sect. [5] we compute 
the running strange quark mass, using the values of Zp obtained from the %SF scheme 
and the tuned bare quark mass obtained using twisted mass Wilson fermions at maximal 
twist. We compare our findings with the strange quark mass obtained using standard SF 
as a renormalization scheme and non-perturbatively improved Wilson fermions. All these 
results are, therefore, a demonstration of the correctness of the continuum limit of the 
%SF renormalization scheme and, consequently, of its applicability in the determination 
of renormalization factors. 

The analysis presented in this paper relies substantially on the results obtained in ref. pp. 
We therefore assume from now on that the reader is familiar with that paper and the 
notation adopted there will be taken over without further notice. Additionally all the 
equations of ref. pQ will be denoted by the equation number prefixed by I as for example 
(1.5.10). 



2 Step-scaling functions: pseudoscalar density 

The SSF, o"o, of a scale-dependent observable, 0(L), describes the behavior of 0(L) under 
changes in the value of the renormalization scale, /x = 1/L, where L denotes the linear 
extent of the finite volume. The reason SSFs are good candidates to perform universality 
tests is that they are finite quantities that depend upon the renormalization scheme and 
the renormalization prescription employed. While at finite values of the cutoff, SSFs are 
regularization-dependent, after the removal of the cutoff they are independent of the 
regulator. 

As a simple example of a SSF, we consider the normalization Zp of the pseudoscalar 
density. The evolution of Zp from a scale L to sL is described by the step-scaling function 
defined by 

Z P (L)a P (s,g 2 (L)) = Z P (sL), (2.1) 

where the g 2 {L) is the renormalized coupling at the scale L. Since the running of the 
coupling is needed to compute the SSF of Zp, or more generally any operator, then the 
SSF of the gauge coupling itself is needed for the computation of any SSF. In this work, 
the quenched setup is chosen and the gauge action is the Wilson gauge action, so there is 
no need to recompute the renormalized gauge coupling and its corresponding SSF. These 
are already known from previous publications [5l6l20] for a very wide range of energies. If 
a different gauge action or dynamical fermions were included, the gauge coupling would 
need to be recomputed with the new formulation. 
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The previous discussion refers to the formal continuum theory. Using a hypercubic lattice 
with spacing a, the chosen renormalization condition for the pseudoscalar density in the 
xSF scheme is 



Z-p(g ,L/a) = c(9,a/L) 



9x{0) 



m=0 



(2.2) 



9p(L/2, 

In this expression, m = indicates that the renormalization condition is imposed at zero 
quark mass. In our case, corresponding to Wilson fermions, this means at the critical value 
of the bare quark mass, mo = m cr , as determined from the tuning in ref. pQ. The factor 
c(9,a/L) is chosen such that Zp takes the correct value at tree-level, Zp(0,L/a) = 1. 
Therefore it is defined as 

c{t), a/L) = 



mo=0 



(2.3) 



The two-point functions entering the definition of Zp have already been discussed in 
ref. pp. They are 

a 3 

gp{x ,9) = g P _{x ,9) = -— 3 
with P 1 denoting the pseudoscalar density and 



(2.4) 



(2.5) 



See ref. pQ for the remaining details of eqs. ( 2.4|2.5 ). In order to determine the renormal- 
ization prescription completely, a value of 6 has to be chosen. In particular, we consider 
here two cases, = (0.5,0.5,0.5) and 6 = (1,0,0). 

Given a fixed value of the renormalization scale, defined through g 2 (L) = u, and a fixed 
value of the lattice spacing, leading hence to a fixed value of L/a, the lattice SSF of the 
pseudoscalar density is given by 



S P (s, u, a/L) 



Zp(g ,sL/a) 



Z P (g ,L/a) 



m=0,g 2 (L)=u 



In the continuum limit, the SSF is finite and takes the value 

<jp(s, u) = lim £p(s, u, a/L) . 

a— >0 



(2.6) 



(2.7) 



From the definition of the SSF, eq. (2.6), it follows that in order to compute the SSF 
at a certain value of the renormalization scale, 1/L, the Z-factor needs to be evaluated 
at both L and sL at the same value of the lattice spacing. In this work we always use 
s = 2. This means that for a fixed value of a (equivalently (3 = 6/g$), simulations have 
to be performed at a certain value of L/a and also at 2L/a. In such computations, the 
values of all parameters, e.g. « cr and zj, are the same at L and 2L for a fixed value of the 
lattice spacing since these parameters only depend on the bare coupling. This is important 
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because it implies that the tuning of the parameters only needs to be performed on the 
smaller lattices. 



We summarize the results for Z-p in tab. A.l at = (0.5, 0.5, 0.5) and = (1, 0, 0). Results 



at three different values of the renormalization scale are presented; the hadronic scale 
L = 1.436 ro, the intermediate scale g 2 = 2.4484 and the perturbative scale g 2 = 0.9944. 
The computation of the SSF is performed only at the intermediate and the perturbative 
scales. We use the results for Zp at the hadronic (matching) scale for the calculation 
of the renormalized strange quark mass as discussed in sect. |5j All results presented in 
the present and following sections have been obtained using the critical values of the 
parameters, k cv and zj, as determined from the tuning condition called (1*) in ref. pp. 



From the values in tab. A.l we have computed the SSF for each lattice spacing at the inter- 
mediate and perturbative scales. The results are presented in tab. A. 2 for 6 = (0.5, 0.5, 0.5) 
and tab. A. 3 for = (1,0,0). In tab. A. 2 we also show the results obtained from the SF 



with improved and standard Wilson fermions [20J. 

For both the %SF and the SF with clover improved Wilson fermions, we have performed 
the continuum limit with a linear fit of the SSF in (a/L) 2 , i.e. using a form 



y = c + ci 



L 



While a linear fit in a/L for the SF with standard Wilson fermions was used, 



y = c + ci 



(2.8) 



(2.9) 



The results of our fits are summarized in tab. A. 4 for 6 = (0.5, 0.5, 0.5) and in tab. A. 5 for 
= (1, 0, 0). In both tables we show the results at the two values of the renormalization 
scale that have been considered, g 2 = 2.4484 and g 2 = 0.9944. For comparison, in tab. A.4 
we also present the continuum limit results for the SF with improved and standard Wilson 
fermions. We have performed our own fits of the values obtained from the SF, since in [20] 
there are no tables with the final continuum limit values, where we could read the data 
from. 

A comparison of the SSF for the different lattice fermions is shown in fig. [T} There we 
plot the results for Ep (cf. tab. A. 2) as a function of a/L for the three formulations. The 



corresponding continuum limit values are also shown (cf. tab. A.4). We can see that the 
slopes, i.e. the values of c\ and c\ in eqs. ( 2.8[2.9 ), are consistent with zero in each case. 
Furthermore the results of the three regularizations agree in the continuum limit at both 
values of the renormalization scale. Moreover, at non-zero lattice spacing the results for 
the three formulations agree at L/a — 16 for the intermediate scale and at Lja — 12, 16 
for the perturbative scale. 



In fig. [2] we show the extrapolation to the continuum limit of £p as determined from the 
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Fig. 1. Lattice SSF of the pseudoscalar density and continuum limit values. Results are shown for 
the xSF with standard Wilson fermions and the SF with improved and standard Wilson fermions, 
at the intermediate and perturbative scales and for 6 = (0.5, 0.5, 0.5). The extrapolations to the 
continuum limit are linear in (a/L) 2 for the xSF and the SF with improved fermions and linear 



in a/L for the SF with standard Wilson fermions. The numerical values are provided in tab. A. 4 
The values from the SF [20] have been slightly displaced to the right and left for the improved 
and the standard Wilson fermions formulations, respectively. 

xSF, now for both values of 6. The results are plotted as a function of (a/L) 2 and the 
corresponding values in the continuum limit, dp, are also shown. The fitting curves are also 
plotted. Note that for the two values of 6 employed in the definition of the renormalization 
prescription, the continuum limit values should not be compared. Different values of give 
rise to different renormalization prescriptions. Consequently, the results are not expected 
to agree in the continuum limit, even though the values appear quite similar. 

A similar plot is shown in fig. |3j where we compare the results of the extrapolation to the 
continuum limit for the xSF and the improved SF. Results are shown at the two values 
of the renormalization scale and for 6 = (0.5,0.5,0.5). Since the same renormalization 
prescription is chosen in both formulations, the continuum limit values should now agree 
between the SF and the x^F. As seen in the plot, this is indeed the case. Note that 
the slopes for the two formulations have similar absolute values but opposite signs. This 
suggests that a constrained continuum limit extrapolation would reduce the uncertainty 
of the final results. 

From the results presented here, we conclude that the continuum limit of the SSF of 
the pseudoscalar density determined from the %SF agrees with that obtained from the 
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Fig. 2. Continuum limit extrapolations of the SSF for the pseudoscalar density. The xSF results 
for both 9\ = 6 = (0.5,0.5,0.5) and 02 = = (1,0,0) are shown at the intermediate and 
perturbative scales. The extrapolations to the continuum limit are linear in (a/L) 2 and plotted 
as dashed lines. The fit results are provided in tab. A. 4 and tab, |A,5 The values in the continuum 



limit are also plotted. The results from 82 have been slightly displaced to the left. 

standard formulation of the SF, both with and without improvement. These results are 
therefore a successful test of the universality of the continuum limit. Moreover, the ap- 
proach to the continuum limit is consistent with the expected leading 0(a 2 ) discretization 
effects in the xSF formulation, indicating that the %SF is indeed compatible with bulk 
automatic 0(a)-improvement. 



3 Step-scaling functions: twist-2 operators 



The previous discussion for the pseudoscalar density can be repeated for any other scale- 
dependent observable. In particular, we now compute the SSF of two lattice realizations 
of the non-singlet twist-two operator whose hadronic matrix elements are related to the 
lowest non-trivial moment of the corresponding unpolarized structure function. 

In Minkowski space the relevant twist-2 gauge-invariant composite operator is 
i n_1 ip(x)j { ^ Df, 2 ■■■ D^} y^( x ) + ' trace terms ' • 
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Fig. 3. Continuum limit extrapolation of the SSF of the pseudoscalar density. Results are shown 
for the xSF with standard Wilson fermions and the SF with improved Wilson fermions, at 
the intermediate and perturbative scales and for = (0.5,0.5,0.5). The extrapolations to the 
continuum limit, plotted as dashed lines, are linear in (a/L) 2 for both cases. The results from 
the fits are presented in tab. A. 4 The values from the SF have been slightly displaced to the 
left. 



The symmetrization in the Lorentz indices, {■ ■ ■}, is required because we deal here with 
unpolarized scattering and the 'trace terms' (terms with g^) are needed in order to 

provide the composite field with a definite spin. The covariant derivative, £) M , is defined 
as the combination 

D,= \ , (3.2) 
with and *D ^ the covariant derivatives acting to the right and left, respectively. 

Since we compute quantities within the %SF scheme and we work in the twisted basis, we 
give here the explicit expressions of the twist-2 operators in Euclidean space and in the 
twisted basis. These can be obtained from the corresponding expressions in the physical 
basis by applying the standard axial rotation eq. (1.2.2) in the continuum theory and then 
directly translating the fields and derivatives to the lattice. This results in 



0;„(x) = x(x) e^ 75T 7{m D„ } - X (x) . (3.3) 
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Depending on the flavor structure, 



x(^)7{m Du} 



cos(a) \ + e ab3 sin(ot) 75 



X (x) (a = 1,2) 



(a = 3) 



(3.4) 



with e a bc the totally anti-symmetric tensor (6123 = 1). In the particular case of maximal 
twist, a = 7r/2, this expression reduces to 



eab3X(x)l{v Dy} 75 \x(?) ( a = 1 ; 2)> 

a 

x(^)7{a« £>*d T^( x ) (a = 3). 



(3.5) 



These composite fields are scale-dependent quantities that need to be renormalized, Or = 
Zq 1 Ob, where Ob is either bare operator from eq. (3.5) and Zq is the corresponding 
renormalization constant. 

The corresponding SSF in the continuum az (s,g 2 (L)) is defined as 

Z (sL) = a Zo (s,g 2 (L))Z (L). (3.6) 

Using a lattice regulator, lattice artifacts must be taken into account and the renormalized 
operator is given a^] 

R (L) = limZo^L/^OBigo). (3.7) 

a— >0 

Eventually, we will compute boundary to bulk correlation functions with the twist-2 op- 
erators inserted in the bulk of the lattice at some space-time point x. The boundary 
interpolating fields at x = that we consider here are, following (21], expressed in the 
physical basis, 



a 6 £C(yh^( z ) 



(3.8) 



Performing a rotation to the twisted basis, with maximal twist angle, such boundary 
interpolating fields take the form 



'e a &3 a 6 Ey,z C(y) 7fc75T Q- C( z ) = 1, 2 ) 



O 



(3.9) 



^ 6 Ev,z C(y)7*T^-C(z) 



(a = 3) . 



2 The definition of the renormalized operator with Z 1 is done to be consistent with the defi- 
nitions used in ref. [7]. 
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In particular, we consider two cases for the gamma matrices at the boundaries, 7^ with 
k = 1, 2. The case k = 1 (k = 2) will be used when computing the correlation function 
of the operator 0£ 4 (0" 2 )- Finally, the correlation functions that we consider here are the 
following, 

£12(2:0, 9) = J2 (x{x) 7{i D 2 } 7s y x(^) C(y) 727sy <3- C(z)> , (3.10a) 

9 11 

£44(^0, 0) = J] 7{o Do} 75 y X(^) C(y) 7i7sy Q- C(z)> • (3.10b) 

Note that we have chosen only the cases where the two flavor matrices, in the bulk and at 
the boundary, are the same and we have picked up only the component r 1 . The reason for 
choosing both matrices to be the same is that, due to symmetries, this is the only combi- 
nation for which the correlation function does not vanish. Amongst the three possibilites, 
r 1,2 ' 3 , each of them should lead to the same value in the continuum limit. However, due 
to our particular setup where flavor symmetry is broken at finite lattice spacing, there is 
a distinction between r 1,2 and r 3 . Choosing r 3 would lead to simpler looking expressions, 
but the appearance of computationally demanding disconnected diagrams leads us to opt 
for t 1 ' 2 . Since these two cases are exactly equivalent, we can select just r 1 . 



We now specify a renormalization prescription for the twist-2 operator within the %SF 
scheme. In particular we impose the renormalization condition 



Z o (g ,L/a) = c(0,a/L) 



9o(L/2,6) 



m=0 



(3.11) 



The factor c(8, a/L) is chosen such that Zq takes the correct value at tree-level, Zq(0, L/a) 
1. Therefore it is defined as 



c(0,a/L) = 



9M 



9o(L/2,9) 



tree 



mo=0 



(3.12) 



In this expression, gi is the two-point function defined in eq. (2.5). The other two-point 
function, g Q , is either g 12 or g 44 from eq. (3.10), corresponding to either 0® 2 or O^, 



respectively. In order to fix the renormalization prescription completely, a value of 6 
has to be chosen. In particular, we consider again the two cases = (0.5, 0.5, 0.5) and 
= (1,0,0). The reason for studying the case with = (1,0,0) is that this is the only 
choice with 0^0 for which calculations with the standard SF have been performed [7]. 
This allows us to compare the continuum limit of the %SF. Although there are no SF 
computations available for the choice 6 = (0.5, 0.5, 0.5), we have also analyzed this setup 
since this value of theta is the usual choice for calculations with the SF formulation. 
Moreover, this additional choice for the parameter 6 allows us to examine the relative 
statistical uncertainties when changing the renormalization prescription through 0. All 
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correlation functions are evaluated at Xq = T/2, where T = L is the time extent of the 
lattice. The scale factor is always set to s = 2. 



Given a fixed value of the renormalization scale, defined through g 2 {L) = u, and a fixed 
value of the lattice spacing a, leading to a fixed value of L/a, the lattice SSF of O, defined 
in the chiral limit, is given as 



s,u,a/L) 



Z o (g ,sL/a) 



m=0,g 2 (L)=u 



Zo(g ,L/a 

In the continuum limit the SSF is finite and takes the value 

R (L) 



ct Zo (s,m) = \imEz Q (s,u,a/L) 



R (sL) 



g 2 (L)=u 



(3.13) 



(3.14) 



with Or(L) the renormalized operator at a given value of the physical scale 1/L computed 
in the xSF renormalization scheme. 

We employ the definitions given above and the chosen renormalization prescription to 
determine the renormalization factors of the operators Of 2 and 0\ A within the %SF scheme 
at finite lattice spacing. Results are presented for several /3-values and at three values of 
the renormalization scale, 1/L. Two values of L correspond to g 2 = 0.9944 and g 2 = 2.4484 
and the third is L = 1.436 tq. These results are given in tab. A. 6 and tab. A. 7 for 0\\ 



and 044, respectively. In both cases we provide results at both 



(0.5,0.5,0.5) and 



= (1,0,0). 



Concerning our results for the Z-factors in the %SF scheme, tab. A. 6 and tab. A. 7 at 
the three values of the renormalization scale and for all values of the lattice spacing that 
we have analyzed, we now make several observations. For either operator 0\ 2 or O44, the 
relative statistical errors in the renormalization constants, AZ0/Z0, are always smaller 
for = (1,0,0) than for 6 = (0.5,0.5,0.5) by nearly a factor of 2. Moreover, at fixed 
values of all parameters, the relative errors in Zq aa are always slightly larger than those 
of Z Ql2 . These results are consistent with the pattern discussed previously in [7] within 
the standard SF setup, where it was shown that for \0\ < 1 the relative statistical errors 
in Zo 12 and Zo 44 increase with decreasing There it was also shown that the statistical 
errors in Zo 4i are slightly larger than those in Zo 12 , which is consistent with what we 
observe from our %SF calculation. 

At the two most perturbative couplings, g 2 = 0.9944 and g 2 = 2.4484, we have de- 



termined the lattice SSFs for both operators. The values are provided in tab. A. 8 for 
6 = (0.5,0.5,0.5) and in tab. A. 9 for 6 = (1,0,0). In tab. A. 9 we have also added 



the values obtained for the lattice SSFs using the SF scheme with standard and non- 
perturbatively improved Wilson fermions. These values were taken from [7] , where results 
are available only at the intermediate coupling, g 2 = 2.4484. 



To take the continuum limit we have performed fits to our results in tab. |A.8|and tab. |A.9 
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Fig. 4. Continuum limit extrapolation of the SSF of the operator 0\2- Only xSF results are 
shown, both for 9\ = 6 = (0.5,0.5,0.5) and 62 = 6 = (1,0,0) and at both the intermediate 
and perturbative scales. The extrapolations to the continuum limit, plotted by dashed lines, are 



linear in (a/L) 2 . The fit results are given in tabs. A. 10 and A. 11 and the values in the continuum 
limit are also plotted. The results for 62 have been plotted slightly displaced to the left. 

linear in (a/L) 2 for the %SF formulation and linear in a/L for the SF with both standard 
and improved Wilson fermions. The fits for the SF calculation are linear in a/L even for 
the formulation with improved fermions because, although the action is improved in this 
setup, the twist-2 operators themselves are not. This is actually the major advantage of 
the xSF as a non-perturbative renormalization scheme: it is not necessary in this setup 
to determine additional counterterms for the operators, since this formulation preserves 
bulk automatic O(a) improvement, up to boundary effects that are expected to be small. 
The results of these fits are presented in tab. |A.1Q and tab. A. 11 



The results for the lattice SSFs from tabs. A. 8 and A. 9 are plotted in figs. [4]{7J where we 
show the continuum limit approach of all SSFs that we have computed. In these figures 
we have also plotted the corresponding values of the SSFs in the continuum limit and the 



fitting curves for the %SF case, as given in tabs. A. 10 and A. 11 



In fig. [4] and fig. [5] we show the continuum limit approach of the SSFs of the operators 0" 2 
and 04 4 , respectively, within the %SF scheme. In each figure, we plot the results for both 
values of 6 and the two scales where we have computed the SSFs. Results are presented 
as a function of (a/L) 2 , since only the %SF computation is considered. We conclude, from 
these figures and the corresponding tables, that the cutoff effects in the %SF SSFs are 
consistent with 0(a 2 ) and are, in fact, small. We can also see that the discretization effects 
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Fig. 5. Continuum limit extrapolation of the SSF of the operator O44. This is the same as fig. [4] 
but now for O44 instead of 0\2- 



are similar for different values of the renormalized coupling and 6. Note that the values 
in the continuum limit for different values of are not required to agree, since different 
values correspond to different renormalization prescriptions. 



In fig. [6] and fig. [7] we compare the results for the SSFs of 0" 2 and 0\ A obtained from 
the three formulations, xSF and SF with standard and improved Wilson fermions. Only 
the computations for the intermediate coupling and for = (1, 0, 0) are plotted. For a 
comparison with the SF, the results are plotted as a function of a/L, although the values 
in the continuum limit for the %SF have been obtained from linear fits in (a/L) 2 . As 
explained earlier, we have performed the continuum extrapolation of the SF results in [7j. 



The results of each of these extrapolations can be found in tab. A. 11 Additionally, the 
continuum limit extrapolations for = (0.5,0.5,0.5) are provided in tab. A. 10 for the 
xSF formulation. There is good agreement, within the statistical uncertainties, between 
the xSF and the improved SF formulations in the continuum limit. 



In summary, we conclude that there is agreement, within the statistical errors, in the 
continuum limit amongst the results from the three formulations. This is another check of 
the universality of the continuum limit, this time through the SSFs of a twist-2 operator. 
Additionally, we observe that the scaling behavior of the SSFs obtained from the %SF is 
consistent with leading 0(a 2 ) discretization effects, which turn out to be rather small. 
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Fig. 6. Continuum limit approach of the SSF of the operator Oyi- Results are shown for the xSF 
with standard Wilson fermions and for the SF with improved and standard Wilson fermions. 
The coupling corresponds to the intermediate scale and = (1,0,0). The continuum limit is 
taken linear in (a/L) 2 for the xSF and linear in a/L for the SF with Wilson and clover improved 



Wilson regularizations. The fit results are provided in tab. ATI] and the continuum limit values 
are also plotted. The data from the SF have been plotted slightly displaced to the right and left, 
respectively, for the improved and unimproved formulations. 

4 Renormalization factors at hadronic scales 



Having established that the %SF does indeed maintain automatic bulk O(a) improvement 
for several SSFs, we now discuss a physical application of these calculations. In this section, 
we determine the renormalization factor at the hadronic (matching) scale for the quark 
mass and for the twist-2 operators discussed in the previous section. Then in the next 
section, we apply this to the determination of the strange quark mass. 

With Wilson twisted mass (Wtm) fermions, the bare vector Ward identity [22] is exactly 
satisfied if one uses a point-split vector current. This implies that Z^ = Zp , where Z^ 
is the bare quark mass renormalization and Zp is the renormalization constant of the 



pseudoscalar density defined in eq. (2.2). In the following, we denote the RGI quark mass 
by M, the bare quark mass is ^ q {go) and Ji q (L) is the renormalized running quark mass 
at the value 1/L of the renormalization scale. The renormalized quark mass is given by 

Jl(L) = Z P \g ,L) N (g ). (4.1) 
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Fig. 7. Continuum limit approach of the SSF of the operator O44. This is the same as fig. [6] but 
now for O44 instead of Oyi- 

The RGI quark mass can be related directly to the bare quark mass by 

M = Z M (g ) N (g ). (4.2) 
The renormalization factor Zu(go) is defined as the product of two terms as follows 

- (m) (dbo) ■ (43) 

The first term, M/]i q (L), is regularization independent but depends on the renormaliza- 
tion scheme as well as on the chosen matching scale 1/L. The second term, Zp 1 (go, L), 
depends on both the renormalization scheme and the regulator. The dependence is such 
that the RGI Z-factor Zu{go) does not depend on the renormalization scheme but only 
on the regularization. All dependence on the matching scale has also disappeared. 

In the discussion above, all equations correspond to the continuum theory. When the 
lattice is used gularization scheme, the correct relation is 

M = Z M (g )ii 9 (g )+O(a n ), (4.4) 

with n = 1 in case of unimproved formulations and n — 2 if improvement is applied. 

The regularization independent part of Zu(go), M/]I q (L), has already been determined 
in [6]. The value is known in the continuum theory and at the matching scale L = 
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1.436 ro- Once the continuum limit is performed, this factor is then universal, i.e. regulator 
independent and we can use it directly in our calculations without the need of a new 
computation since both the SF and the xSF are equivalent formulations in the continuum 
theory. The value obtained in [6] is 



M/jL q (L) = 1.157(15) at L = 1.436 r 
with a relative error of 1.3%, which is sufficient for our purposes. 



(4.5) 



This means that we are left with the computation of two quantities. One is the regu- 
larization dependent part of the total renormalization factor, Zp 1 (go, L), which must be 



computed at the matching scale L of eq. (4.5) for several values of /3 and within the 
xSF scheme. The other quantity is the bare quark mass, fi q (go), which also has to be 
determined for a range of bare couplings and is discussed shortly in sect. [5] 

The determination of Zp(g , L/a) from the x^F, at a certain value of the renormalization 
scale and for several values of the lattice spacing, has already been explained in sect. [2] 
and the results are given in tab. A.l| Amongst the cases presented in tab. A.l we can 
restrict our focus to the results corresponding to the matching scale, L = 1.436 tq and 
6 = (0.5,0.5,0.5). For this choice of the parameters, we have computed Zp(g , L/a) at 



several values of the lattice spacing in the range 6.0 < (3 < 6.5, which we recall in tab. B.l 



With these results, we determine a smooth parameterization of the dependence of Zp on 
/3 for the range covered in our calculation. We use a polynomial fit to parameterize our 
results, 



Z P (g ,L/a) L=1A36ro = ^2zJ(f3 - 6.0) 4 

i=0 

f3 = Q/g 2 o, 6.0 < /3 < 6.5. 



(4.6) 



The fitted coefficients are provided in tab. B.2 



This paramterization of Zp, combined with eqs. (4.3) and (4.5), provides a parameteriza- 
tion of Zm itself. It has the form 



i=0 

(3 = G/g 2 , 6.0</3<6.5. 
with the coefficients presented in tab. 



(4.7) 



B.2 



The uncertainty of M/~p q (L) is independent 
of (5 and hence can be accounted for after the extrapolation to the continuum limit has 
been carried out. 



From eq. (4.6) and eq. (4.7), it is now possible to compute Zp and Zm at any value of (3 



within the range 6.0 < (3 < 6.5. This is the range of (3 where large volume calculations 
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have been performed, namely f3 = 6.00, 6.10, 6.20, 6.45. For these computations, a number 
of bare quark masses were used, which we can use to determine the renormalized quark 
mass from our knowledge of Z P and Z M . The relevant values of Z P and Z M at the chosen 



values of (3 are summarized in tab. B.4 



As we just did for the pseudoscalar density, we now determine the RGI Z-factors of the 
operators Of 2 an d 0% 4 using the xSF formulation. These factors then relate the bare and 
the RGI matrix elements of the corresponding operator. 

We first study the dependence of the Z-factors on (3 at the matching scale L = 1.436 ro 
and determine a curve describing this dependence. Performing a fit of the values for the 



Z-factors (found in tabs. A. 6 and A. 7) of the form 



Z o (g , L/a) L=1A36ro = £ z REN (/3 - 6.0) 4 , 

i=o l 4 - 8 J 
(3 = 6/g 2 , 6.0</3<6.5, 



we obtain the fitting coefficients z REN given in tab. B.5 for the %SF and the SF with 



standard and improved Wilson fermions. In eq. (4.8), 'REN' stands for the particular setup 
chosen: xSF, SF with standard Wilson fermions or SF with improved Wilson fermions. We 
have computed the Z-factors for 0" 2 and 0£ 4 at 6 = (1,0,0) for all formulations and at 
= (0.5, 0.5, 0.5) only for the x^F. These results are obtained from fits performed in this 
work for the three formulations, which for the SF are in agreement with the final results 
previously presented in [23]. We show the Z-factors together with the fitting curves, at 
6 = (1, 0, 0) and for the three formulations, in figs. M and M for Zq 12 and Zq 44 , respectively. 



As for the case of the pseudoscalar density, it is possible to determine the RGI renor- 
malization constants for the twist-2 operators from the knowledge of Z REN (go, L/a) at a 
given value of the renormalization scale, 1/L. For the twist-2 operators one defines, at the 
same value of the renormalization scale, the ultraviolet (UV) invariant SSF, <tJ^ i ^ n (L). 
This SSF is discussed thoroughly in ref. [24], which we assume some familiarity with. 



The UV invariant SSF of a certain operator is independent of the particular regularization 
but it depends on the renormalization scheme and the matching scale. In particular it is 
defined by 

RGI = Pa(L) , (4.9) 



with Or(L) the renormalized operator at scale 1/L, defined in eq. (3.7), and the corre- 
sponding RGI operator RGI . We note that the UV invariant SSF is analogous to the 
factor \M /~p q (L) \ that we have used for the RGI quark mass. 

The RGI renormalization factor is scale and scheme independent but it depends on the 
particular regularization. It relates any bare matrix element of the bare operator, Ob (go); 
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Fig. 8. Results for Zo 12 (go, L/a) at the scale L = 1.436ro, = (1,0,0) and for several values 



of (3. Results are shown for the xSF with standard Wilson fermions (cf. tab. A. 6) and for the 
SF with standard and improved Wilson fermions, as taken from [23] . The fitting curves are also 
plotted (cf. eq. ([48]) and tab. pO])). 

with the corresponding RGI matrix element and it is defined as follows, 



7 rgi 
J o 



(9o) 



7 REN 



(9o, L/a) 



a 



UV.REN 
INV.O 



(L) 



(4.10) 



In [7J, the value of the UV invariant SSF was given for the operators 0\ 2 and 0^4 a ^ the 
scale L = 1.436 r and for 6 = (1, 0, 0). The values given there are 



a 



UV,SF 
INV,Oi2 



0.242 (8) 



a 



UV,SF 
INV,0 4 4 



0.221 (9) . 



(4.11) 



Substituting these values into eq. (4.10) and using the Z-factors, at the matching scale 
L = 1.436 ro, in tabs. A. 6 and A. 7 the RGI Z-factors of the operators 0\ 2 and are 
determined and given in tab. B.6 Results are provided only for 6 = (1,0,0) because 
the UV invariant SSFs are only known for that case. In the determination of the RGI 
Z-factors, the uncertainty in the UV invariant SSF is not taken into account. This is a 
quantity in the continuum, and therefore its uncertainty is only considered at the end of 
all calculations, after the continuum limit has been performed. Its error is then added in 
quadrature to the final uncertainty in the continuum limit. 



We smoothly parameterize Zq (go) as a function of (5 by fitting the values in tab. B.6 to 
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Fig. 9. Results for Zo 44 (go, L/a) at the scale L = 1.436 ro, 6 
p. This is the same as fig. [8] but now for O44 instead of 0\2- 

the following functional form 



(1,0,0) and for several values of 



^o GI M = E^ GI (/ 3 - 6 -°r 

i=0 

(3 = 6/g 2 , 6.0</3<6.5. 
The fitted coefficients can be found in tab. IB. 71 



(4.12) 



Using the parameterizations in eqs. (4.8) and (4.12), both the Z-factors and the RGI Z- 



factors can be determined at any value of f3 within the range 6.0 < (3 < 6.5. In tab. B.8 



we provide results for the particular values of f3 for which bare matrix elements have been 
evaluated in large volume calculations [25] . 



The results presented in tab. B.8 correspond to Zq 12 , Zq^, Zq 44 and Zq^ at 6 = (1, 0, 0). 



We also give there the corresponding results for the SF formulation with improved and 
standard Wilson fermions. Note that these values should not be compared across the 
three formulations. They depend on the regularization. Only a comparison of renormalized 
matrix elements in the continuum limit would make sense. 



Nevertheless, this calculation of the renormalization constants demonstrates that the %SF 
can indeed be used in pratice to non-perturbatively renormalize challenging operators, 
such as the twist-two operators considered in this work. 
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5 Strange quark mass 



In this section we compute the RGI strange quark mass, M g , and the running strange quark 
mass in the MS-scheme at 2 GeV in quenched QCD. We use the %SF renormalization 
scheme with the setup discussed in the previous sections together with the bare quark 
masses from large volume calculations with twisted mass fermions at maximal twist. 

The purpose of this computation is to perform another check of the xSF formulation. In 
practice, we compute the quantity r (M s + M), where M = (M u + Md)/2 is the average 
light quark mass and r is the Sommer parameter [26J. We then take the continuum limit. 
The resulting continuum limit value, obtained from the xSF, is compared to that obtained 
from the standard SF with improved Wilson fermions [21] ■ We find that the two results 
agree, which is another check of universality in the continuum limit. 

Moreover, r$ (M s + M) is expected to scale towards the continuum limit with leading 
0(a 2 ) discretization errors, up to possible boundary effects. In fact, we will show that the 
scaling behavior is consistent with leading 0(a 2 ) discretization effects. This represents 
another test of bulk automatic 0(a) improvement and, moreover, it provides another 
indirect indication that the boundary effects coming from d s (see ref. pQ) are negligible, 
even at the large values of go considered in this section. 

To determine the strange quark mass, we follow the strategy of ref. [21] and determine 
a reference bare quark mass /i re f defined by 2/i re f = fi s + fi, where fi = (fi u + /Ud)/2. The 
reference quark mass is then chosen such that the physical value of the kaon meson mass 
is reproduced. 

To start, we determine the bare reference quark mass in lattice units, a/z re f (go), using the 
results of the large volume calculations in [28]. The pseudoscalar mass, m PS , was computed 
there in lattice units using twisted mass Wilson fermions at maximal twist. We focus on 
the values of /3 in [2E] that overlap with the range covered in this work 6.0 < (3 < 6.5. 

The pseudoscalar mass range covered by these data set is 270 MeV < m PS < 1180 MeV. 
Within this range, we may interpolate in the bare quark mass a/j, q at the experimental 
value of the kaon mass, m^. At each value of 0, we perform a quadratic interpolation in 
the bare quark mass. We have cross checked that a linear interpolation with the 3 data 
points closest to the interpolation point gives consistent results. The lattice spacing in 
physical units is obtained using the dependence of r /a from ref. [23] • The final results 
for /iref^o, together with the corresponding values of r /a and a/i re f, are summarized in 
tab. ED 

Using the results for Zm((?o) obtained in sec. [5] and the reference quark mass just dis- 
cusssed, we can determine the RGI strange quark mass. The resulting values at finite 
lattice spacing for the RGI reference quark mass are summarized in tab. [2] and plotted in 
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p 


ro/a 


a/i r ef 


/^rcf^O 


Kcr° n definition 


6.00 


5.368 (22) 


0.01450(59) 


0.0778(32) 


6.10 


6.324(28) 


0.01216(40) 


0.0769(26) 


6.20 


7.360 (35) 


0.01030(34) 


0.0758(25) 


6.45 


10.458 (58) 






K cr CAC definition 


6.00 


5.368(22) 


0.01443(51) 


0.0775 (28) 


6.20 


7.360 (35) 


0.01029(27) 


0.0757(20) 



Table 1 

Hreffo at the value of the kaon mass and for particular values of /3. Results are shown for both 
the pion and the PCAC definitions of the critical mass, as discussed in ref. |28j . 



p 


Zu 


r (M s + M) 


Kcr° n definition 


6.00 


2.1444 (55) 


0.334(14) 


6.10 


2.1733 (33) 


0.334(11) 


6.20 


2.1957 (42) 


0.333(11) 


6.45 


2.2236 (70) 




K cr CAC definition 


6.00 


2.1444 (55) 


0.332(12) 


6.20 


2.1957 (42) 


0.3324(88) 



Table 2 

ro (M s + M) = Z]vi(2// re f ro) at several values of (3. Results are shown for the pion and the PCAC 
definitions of the critical mass. 



fig. 10 



We plot the results for r (M s + M) obtained using the two methods for determining K cr 
in ref. The subtleties related to these 2 different choices are discussed in ref. [30] 
and refs. therein. For this work, these two definitions simply correspond to two slighlty 
different discretizations of the twisted mass action inducing slightly different 0(a 2 ) cutoff 
effects in the physical quantities. 



In fig. 10 we also show a continuum extrapolation, linear in (a/ro) 2 , and the resulting con- 
tinuum limit values. The results obtained from the SF with improved Wilson fermions [27] 
are also plotted for comparison. The final values in the continuum limit are given in 
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Fig. 10. Scaling of vq(M s + M) at the physical value of the kaon mass, dik = 495 MeV. The 
extrapolations to the continuum limit are performed with linear fits in (a/ro) 2 . The values in 
the continuum limit are also plotted. Results are shown for the xSF with standard Wilson 
fermions, for the two definitions of the critical mass, and also for the SF with improved Wilson 
fermions [27] • The results for the %SF have been plotted slightly displaced to the right and left, 
respectively, for the PCAC and pion definitions of the critical mass. 



eqs. (5.1) and (5.2) for the pion and PCAC definitions of the critical mass, 



pion : r (M s + M) xSF = 0.332 (28) , 
PCAC : r (M s + M) xSF = 0.333 (23) . 



(5.1) 
(5.2) 



The errors include the uncertainty of Mj[i q {V) added in quadrature. These values are 
consistent with that obtained in [27] 



r (M s + M) SF = 0.362 (12 



(5.3) 

As can be seen from these results and those in fig. [TUJ the SF values have relative errors 
that are about two times smaller than those obtained in our calculation using the %SF. 
This difference is due to the size of the statistical uncertainties in the bare pseudoscalar 
masses obtained from the corresponding large volume calculations and not due to an 
inherit difference in the accuracies that can be obtained with %SF or SF computations. 
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Our xSF values rely on the results of [28] whereas the SF values use those in [27], which 
are about twice as accurate as the values in [2H] • This accounts for the larger uncertainty 
of the xSF results for the strange quark mass. 

We can conclude that the values of the RGI reference quark mass, and therefore the RGI 
strange quark mass itself, determined using the SF and the %SF agree in the continuum 
limit. This is another test of the universality, this time at a rather large value of the 
physical volume. In particular, these results demonstrate that the x^F, like the SF, could 
also be a valuable tool for the computation of the renoramalized quark masses. 

Even if not necessary for testing universality, we compute, for completeness, the values 
of both the RGI strange quark mass in physical units and the running strange quark 
mass in the MS-scheme. As discussed in [27], chiral perturbation theory allows for a 
precise determination of ratios of the masses of the three lightest quarks |31|l32U33j . Such 
determinations lead to 

M u /M d = 0.553 ± 0.043 , M B /M = 24.4 ± 1.5 . (5.4) 



Assuming these relations together with eqs. (5.1) and (5.2), we can determine the value 



of the RGI strange quark mass. The final results in units of ro are 

pion : r M* SF = 0.319 (27) , (5.5) 
PCAC : r Mf F = 0.320 (22) , (5.6) 

for the pion and PCAC definitions of k ct , respectively. Repeating this analysis for the 
results in [27] we obtain, 

r M s SF = 0.348(12). (5.7) 

The value of the RGI strange quark mass can be now given in physical units, 

pion : M* SF = 126 (11) MeV , (5.8) 

PCAC : M* SF = 126 (9) MeV , (5.9) 

and for the SF, 

M S SF = 137 (5) MeV . (5.10) 

Using the conversion factor between the RGI mass and the running mass in the MS- 
scheme, the running strange quark mass in the MS-scheme can be determined. At a value 
of the energy scale of 2 GeV and up to 4 loops, the flavor- independent conversion factor 
is [27] 

m^(2 GeV) /M = 0.72076 . (5.11) 
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As a result, the strange quark mass at 2 GeV with 4-loop running in the MS-scheme is 



where we quote as our best result the strange quark mass obtained from the PCAC 
definition of k ct . 



6 Conclusions 

Presently to renormalize Wilson twisted mass fermions preserving the property of au- 
tomatic 0(a) improvement, infinite volume renormalization schemes such as the RI- 
MOM [34.35J or the X-space schemes [36|37] are used. However, only finite volume schemes 
such as the SF [2.38J and the xSF [15] solve the problem of covering a large range of scales. 
Furthermore, the xSF scheme, when used to renormalize bare matrix elements computed 
with maximally twisted mass fermions, is also compatible with automatic 0(a) improve- 
ment. 

In this work we have made a detailed study of several applications of the %SF scheme, 
with quenched Wilson fermions. Using our results for the non-perturbative tuning of the 
xSF [T], we have performed a number of scaling studies of the %SF. We have analyzed 
step-scaling functions of the pseudoscalar density and of two discretizations of a twist- 
2 operator at a perturbative and an intermediate value of the renormalized coupling. 
All our results agree in the continuum limit with those obtained from the standard SF 
scheme. Additionally, our results are consistent with scaling violations of only 0(a 2 ), thus 
demonstrating the bulk 0(a) improvement of the %SF scheme. 

We remark that for the twist-2 operators this is an important result because to improve 
such operators in the standard SF scheme would require a non-perturbative computation 
of additional improvement coefficients. The automatic 0(a) improvement found here with 
the example of the twist-2 operators can be taken over to other observables. It demon- 
strates that automatic 0(a) improvement is at work and that with the %SF scheme, in 
combination with maximally twisted mass fermions, the somewhat demanding computa- 
tion of operator specific improvement coefficients can be avoided. 

Additionally, we have computed the continuum limit of the renormalized strange quark 
mass within the x^F scheme. In this well, the result in the continuum limit 

is consistent with previous results obtained with non-perturbatively improved Wilson 
fermions and the standard SF scheme, thus demonstrating that the %SF scheme works as 
well as the standard SF scheme even for the more commonly computed quantities, such 
as the running of the quark masses. 



XSF: /Is ( 2 GeV ) = 91 ( 6 ) MeV > 
SF: mf (2 GeV) = 99 (4) MeV , 



(5.12) 
(5.13) 
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Therefore, we believe that the scheme is both a practical and theoretically well 
defined framework that can be used as a non-perturbative renormalization scheme for 
large volume calculations of bare operators. 
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A Tables of numerical results for the step-scaling functions 





= (0.5,0.5,0.5) 


e = (1,0,0) 


L/a 


P 


Zp(go,L/a) 


Z P (g ,2L/a) 


Zp(go,L/a) 


Z P (g ,2L/a) 


Hadronic scale: L = 1.436 ro 


8 


6.0219 


0.5385(12) 




0.5432 (12) 




10 


6.1628 


0.5264(12) 




0.5310 (12) 




12 


6.2885 


0.5272(16) 




0.5321 (17) 




16 


6.4956 


0.5187(22) 




0.5245 (21) 








Intermediate scale: g 2 = 


2.4484 




8 


7.0197 


0.68509 (95) 


0.6199(14) 


0.68850 (93) 


0.6241 (13) 


12 


7.3551 


0.6735 (13) 


0.6082(19) 


0.6788 (12) 


0.6142(21) 


16 


7.6101 


0.6672 (16) 


0.5991 (22) 


0.6737(16) 


0.6015(22) 






Perturbative scale: g 2 = 


0.9944 




8 


10.3000 


0.82689 (56) 


0.80129 (84) 


0.83007(58) 


0.80358 (84) 


12 


10.6086 


0.81651 (88) 


0.78549 (84) 


0.81924 (82) 


0.79008 (80) 


16 


10.8910 


0.8110 (10) 


0.7820 (14) 


0.8136(11) 


0.7802 (15) 



Table A.l: Rcnormalization factors of the pseudoscalar density, Zp, at 6 = (0.5,0.5,0.5) and 6 = (1,0,0). Results are 
shown for the xSF with standard Wilson fermions at three values of the renormalization scale and for several values of the 
lattice spacing. 



£p(2,«,a/L) 



L/a 


X SF 


SF (Clover) 


SF (Wilson) 


Intermediate scale: g 2 = 2.4484 


8 


0.9048 (23) 


0.8945 (23) 


0.8993 (20) 


12 


0.9030 (33) 


0.8908 (23) 


0.8924 (30) 


16 


0.8980 (39) 


0.8998 (25) 


0.9036 (32) 


Perturbative scale: g 2 = 0.9944 


8 


0.9690(12) 


0.9633 (14) 


0.9641 (12) 


12 


0.9620(15) 


0.9599 (19) 


0.9632 (17) 


16 


0.9643 (22) 


0.9622 (20) 


0.9652 (22) 



Table A. 2: SSF of the pseudoscalar density at finite lat- 
tice spacing, Sp (2, u, a/L), for 6 = (0.5,0.5,0.5). Results are 
shown for the xSF with standard Wilson fermions and also for 
the SF with improved and standard Wilson fermions 1201 at 
two values of the renormalization scale and for several values 
of the lattice spacing. 



S P (2,u, a/L) 


L/a 


X SF 


Intermediate scale: g 2 = 2.4484 


8 


0.9065 (23) 


12 


0.9048 (35) 


16 


0.8929 (39) 


Perturbative scale: g 2 = 0.9944 


8 


0.9681 (12) 


12 


0.9644 (14) 


16 


0.9590 (23) 



Table A. 3: SSF of the pseudoscalar density at finite lattice 
spacing, £p(2, u, a/L), for = (1,0,0). Results are shown 
for the xSF with standard Wilson fermions at two values of 
the renormalization scale and for several values of the lattice 
spacing. 
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X SF 


SF (Clover) 


SF (Wilson) 






XSF 


Intermediate scale: g 2 = 2.4484 




Intermediate scale: g 2 = 2.4484 


op (2, u) 


0.8981 (41) 


0.8968 (28) 


0.8993 (58) 




Op (2, u) 


0.8942 (42) 


slope 


0.44(34) 


-0.22(27) 


-0.007(56) 




slope 


0.82(34) 


X 2 /dof 


0.5349 


6.4083 


6.8156 




X 2 /dof 


3.3471 


Perturbative scale: g 2 = 0.9944 




Perturbative scale: g 2 = 0.9944 


Op (2, u) 


0.9595(21) 


0.9602 (22) 


0.9644 (37) 




crp (2, u) 


0.9591 (21) 


slope 


0.59(18) 


0.19(19) 


-0.004 (35) 




slope 


0.59(17) 


X 2 /dof 


2.4748 


1.1258 


0.5129 




X 2 /dof 


1.8644 



Table A. 4: Continuum limit of the SSF of the pseudoscalar 
density. Results are shown for the xSF with standard Wil- 
son fermions and also for the SF with improved and standard 
Wilson fermions at two values of the renormalization scale 
and for = (0.5,0.5,0.5). These results correspond to linear 
fits of the values in tab. |A,2"] The fit is linear in a/L for the 
SF(Wilson) formulation while it is linear in (a/L) 2 for the xSF 
and SF(Clover) formulations. 



Table A. 5: Continuum limit of the SSF of the pseudoscalar 
density. Results are shown for the xSF with standard Wilson 
fermions at two values of the renormalization scale and for 
= (1,0,0). These results correspond to linear fits in (a/L) 2 
of the values in tab. |A.3| 





= (0.5,0.5,0.5) 


= (1,0,0) 


L/a 


P 


z Oi 2 i90,L/a) 


Zoi 2 (90,2L/a) 


z O 12 {90,L/a) 


Z Ol2 (g ,2L/a) 


Hadronic scale: L = 1.436 rg 


8 


6.0219 


0.395(12) 




0.3746 (59) 




10 


6.1628 


0.374(13) 




0.3509(59) 




12 


6.2885 


0.348(15) 




0.3547(77) 




16 


6.4956 


0.353(21) 




0.341 (11) 








Intermediate scale: g 2 = 


2.4484 




8 


7.0197 


0.6077 (80) 


0.482(13) 


0.5675 (41) 


0.4498 (64) 


12 


7.3551 


0.613(11) 


0.495(16) 


0.5634(58) 


0.4401 (84) 


16 


7.6101 


0.611 (14) 


0.460(16) 


0.5587 (70) 


0.4383 (84) 






Perturbative scale: g 2 = 


0.9944 




8 


10.3000 


0.7989 (44) 


0.7570 (68) 


0.7717(25) 


0.7287(36) 


12 


10.6086 


0.7800 (66) 


0.7597 (62) 


0.7530 (35) 


0.7295 (34) 


16 


10.8910 


0.7762(83) 


0.730(11) 


0.7511 (45) 


0.7161 (68) 



Table A. 6: Renormalization factors Zq 12 for = (0.5,0.5,0.5) and = (1,0,0). Results are shown for the xSF with 
standard Wilson fermions at three values of the renormalization scale and for several values of the lattice spacing. 
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= (0.5,0.5,0.5) 


= (1,0,0) 


L/a 


P 


Zq 44 (go, L/a) 


Z Oii (90, 2L/a) 


Zo 44 (go,L/a) 


Z Oii (g ,2L/a) 


Hadronic scale: L = 1.436 rg 


8 


6.0219 


0.319(10) 




0.3416 (62) 




10 


6.1628 


0.307(10) 




0.3305 (62) 




12 


6.2885 


0.280(14) 




0.3217(86) 




16 


6.4956 


0.261 (19) 




0.297(12) 








Intermediate scale: g 2 = 


2.4484 




8 


7.0197 


0.5174(75) 


0.388(12) 


0.5382 (44) 


0.4203 (71) 


12 


7.3551 


0.532(11) 


0.404(16) 


0.5340 (64) 


0.4189 (93) 


16 


7.6101 


0.521 (14) 


0.405(17) 


0.5236 (82) 


0.4417(86) 






Perturbative scale: g 2 = 


0.9944 




8 


10.3000 


0.7369 (46) 


0.6781 (71) 


0.7529 (26) 


0.7044 (39) 


12 


10.6086 


0.7145 (64) 


0.6882 (68) 


0.7334 (37) 


0.7068 (39) 


16 


10.8910 


0.7114(88) 


0.686(11) 


0.7301 (51) 


0.7135(65) 



Table A. 7: Renormalization factors Zq 44 for = (0.5,0.5,0.5) and 8 = (1,0,0). Results are shown for the xSF with 
standard Wilson fermions at three values of the renormalization scale and for several values of the lattice spacing. 



XSF 


L/a 


S 0l2 (2,M,a/L) 


Eo 44 (2,u,o/L) 


Intermediate scale: 


g 2 = 2.4484 


8 


0.793 (24) 


0.751 (26) 


12 


0.809 (31) 


0.761 (34) 


16 


0.752(31) 


0.777 (39) 


Perturbative scale: 


g 2 = 0.9944 


8 


0.948 (10) 


0.920(11) 


12 


0.974(11) 


0.963(13) 


16 


0.940 (18) 


0.964 (20) 



Table A. 8: SSF of O12 and O44 at finite lattice spacing, Sq 12 and Sq 44 . Results are shown for the xSF with standard 
Wilson fermions at two values of the renormalization scale, for several values of the lattice spacing and for = (0.5, 0.5, 0.5). 
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So 12 (2,u,o/i) 


So 44 (2,n,a/L) 


L/a 


XSF 


SF (Clover) 


SF (Wilson) 


XSF 


SF (Clover) 


SF (Wilson) 






Intermediate scale: 


g 2 = 2.4484 






8 


0.793(13) 


0.8223 (77) 


0.8811 (85) 


0.781 (15) 


0.7885 (91) 


0.7935(119) 


12 


0.781 (17) 


0.8053 (77) 


0.8589 (136) 


0.784 (20) 


0.7921 (94) 


0.7942(186) 


16 


0.785(18) 


0.8116(107) 


0.8519 (85) 


0.844 (21) 


0.8036(127) 


0.7823(115) 






Perturbative scale: 


g 2 = 0.9944 






8 


0.9443 (56) 






0.9355 (62) 






12 


0.9688 (64) 






0.9637(72) 






16 


0.953(11) 






0.977(11) 







Table A. 9: SSF of O12 and O44 at finite lattice spacing, S012 ano ^ ^044- Results are shown for the %SF with standard 
Wilson fermions and also for the SF [7] with improved and standard Wilson fermions at two values of the renormalization 
scale, for several values of the lattice spacing and for 6 = (1,0,0). 



X SF 




O12 


O44 


Intermediate scale: g 2 = 2.4484 


cro (2, u) 
slope 
X 2 /d°f 


0.766(35) 
2(3) 
1.4081 


0.779 (42) 
-2(4) 
0.0421 


Perturbative scale: g 2 = 0.9944 


cro (2, u) 
slope 
X 2 /d°f 


0.970(16) 
-1(1) 
3.3364 


0.989(19) 
-4(2) 
0.2731 



Table A. 10: Continuum limit of the SSF, (ro 12 and (ro 44 , of the operators O12 and O44. Results are shown for the xSF 
with standard Wilson fermions at two values of the renormalization scale and for = (0.5, 0.5, 0.5). These results correspond 

The fits are linear in (a/L) 2 . 



to linear fits of the values in tab. 
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O12 


O44 




XSF 


SF (Clover) 


SF (Wil) 


XSF 


SF (Clover) 


SF (Wil) 






Intermediate scale: g 2 


= 2.4484 








0.778 (20) 


0.790(19) 


0.822(18) 


0.837 (23) 


0.812(23) 


0.774 (25) 


slope 


1(2) 


0.25(19) 


0.47(19) 


-4(2) 


-0.19(23) 


0.17(26) 


X 2 /dof 


0.0772 


0.8338 


0.0332 


2.9018 


0.2471 


0.1586 






Perturbative scale: g 2 


= 0.9944 






co(2,«) 


0.9752 (97) 






0.988(10) 






slope 


-1.89(82) 






-3.40(88) 






X 2 /dof 


2.9783 






0.0534 







Table A. 11: Continuum limit of the SSF, <jq 12 and o"o 44 , of the operators O12 and O44. Results are shown for the xSF 
with standard Wilson fermions and also for the SF with improved and standard Wilson fermions at two values of the 
renormalization scale and for 6 = (1,0,0). These results correspond to linear fits of the values in tab. |A.9| The fits are 
linear in (a/L) 2 for the xSF formulation while they are linear in a/L for the SF. 
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B Tables of numerical results for the renormalization factors 



L/a 


P 


Z P 


8 


6.0219 


0.5385(12) 


10 


6.1628 


0.5264(12) 


12 


6.2885 


0.5272(16) 


16 


6.4956 


0.5187(22) 



Table B.l: Z P (g ,L/a) at L = 1.436 r and for = 
(0.5,0.5,0.5). 



i 


«f 


z i 





0.5394(14) 


2.1444(55) 


1 


-0.077(15) 


0.321 (60) 


2 


0.078(30) 


-0.32(12) 



Table B.2: Coefficients of the beta dependence of 
Zp(gQ,L/a) at the matching scale L = 1.436 ro (cf. Eq. ( |4.6[ l) 
and Z M (go) (cf. Eq. p~7)). Results are shown for the xSF with 
standard Wilson fermions at = (0.5,0.5,0.5). 



L/a 


P 




Zm 




P 


z F 


Zm 


8 


6.0219 


1.8569 (41) 


2.1484 (47) 




6.00 


0.5394(14) 


2.1444(55) 


10 


6.1628 


1.8995 (42) 


2.1977(49) 




6.10 


0.53240(82) 


2.1733(33) 


12 


6.2885 


1.8968 (59) 


2.1946(68) 




6.20 


0.5270(10) 


2.1957(42) 


16 


6.4956 


1.9279 (80) 


2.2306(93) 




6.45 


0.5203 (17) 


2.2236(70) 



Table B.3: Z~ 1 (g , L/a) at L = 1.436 r andZ M (go), for = Table B.4: Z P (g ,L/a) at L = 1.436 ro and Z M (g ), for = 

(0.5, 0.5, 0.5). Results are shown for the xSF formulation and (0.5, 0.5, 0.5). Results are presented for the xSF formulati on at 

for all the values of P where simulations have been performed. several value s of P, as determined from the curves in eq. l |4.6| l 

and eq. \i.7\ . 





Zi(0 12 ) 


Zi(OAi) 


i 


XSF 


SF (Clover) 


SF (Wilson) 


X SF 


SF (Clover) 


SF (Wilson) 


= (0.5,0.5,0.5) 





0.402 (14) 






0.323(12) 






1 


-0.25(15) 






-0.12(12) 






2 


0.30(29) 






-0.02(25) 






e = (1,0,0) 





0.3761 (69) 


0.3410(31) 


0.3659 (35) 


0.3426 (73) 


0.3450 (37) 


0.3197(44) 


i 


-0.151 (72) 


-0.077(31) 


-0.102(35) 


-0.059 (77) 


-0.180(37) 


-0.117(44) 


2 


0.18(15) 


0.061 (62) 


0.047(70) 


-0.07(16) 


0.196(72) 


0.105(89) 



Table B.5: Coefficients of the beta-dependence of Zq 12 and 2q 44 at the matching scale L = 1.436 ro- Results arc shown 
for the xSF with standard Wilson fermions at = (1,0,0) and = (0.5,0.5,0.5) and also for the SF 1231 with improved 
and standard Wilson fermions at = (1, 0, 0). 
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L/a 


P 


XSF 


SF (Clo) 


SF (Wil) 


XSF 


SF (Clo) 


SF (Wil) 


8 


6.0219 


1.548 (24) 


1.414(13) 


1.518(14) 


1.546(28) 


1.562 (17) 


1.453 (20) 


10 


6.1628 


1.450(24) 


1.347(13) 


1.433(14) 


1.495(28) 


1.450(16) 


1.355(19) 


12 


6.2885 


1.466 (32) 


1.358(12) 


1.431 (14) 


1.456(39) 


1.417(16) 


1.361 (19) 


16 


6.4956 


1.409 (45) 


1.309(15) 


1.348(18) 


1.344(54) 


1.371 (19) 


1.295(26) 



Table B.6: RGI renormalization factors, Z^(g ) and Zg GI (g ) for = (1,0,0). Results are shown for the xSF with 
standard Wilson fermions and the SF with improved and standard Wilson fermions and for several values of the lattice 
spacing. We have determined in this work the RGI Z-factors for the SF from the Z-factors given in |23| . 





-zf GI (0 12 ) 


*f G, (0 4 4) 


i 


X SF 


SF (Clover) 


SF (Wilson) 


X SF 


SF (Clover) 


SF (Wilson) 





1.554(29) 


1.409(13) 


1.512(14) 


1.550 (33) 


1.561 (17) 


1.447(20) 


1 


-0.62(30) 


-0.32(13) 


-0.42(14) 


-0.27(35) 


-0.81 (17) 


-0.53(20) 


2 


0.74(62) 


0.25(26) 


0.19(29) 


-0.32(72) 


0.89 (33) 


0.48 (40) 



Table B.7: Coefficients of the beta-dependence of Zq GI (9o) and Zq GI (<;o). Results are shown for the xSF with standard 
Wilson fermions and the SF with improved and standard Wilson fermions at 6 = (1, 0, 0). 



p 


z Ql2 




Z Ou 






X SF 


6.00 


0.3761 (69) 


1.554(29) 


0.3426 (73) 


1.550 (33) 


6.10 


0.3627(40) 


1.499(17) 


0.3361 (43) 


1.521 (19) 


6.20 


0.3529 (49) 


1.458(20) 


0.3283 (53) 


1.486(24) 


6.45 


0.3436 (82) 


1.420(34) 


0.3031 (90) 


1.371 (41) 




SF (Clover) 


6.00 


0.3410(31) 


1.409(13) 


0.3450 (37) 


1.561 (17) 


6.10 


0.3338(20) 


1.3793 (83) 


0.3290 (23) 


1.489(10) 


6.20 


0.3279(23) 


1.3550 (95) 


0.3168(27) 


1.433(12) 


6.45 


0.3186 (30) 


1.317(12) 


0.3035 (35) 


1.373(16) 




SF (Wilson) 


6.00 


0.3659 (35) 


1.512(14) 


0.3197 (44) 


1.447(20) 


6.10 


0.3562 (22) 


1.4719(91) 


0.3091 (28) 


1.399(13) 


6.20 


0.3474 (25) 


1.436(10) 


0.3005 (32) 


1.360(14) 


6.45 


0.3296 (36) 


1.362(15) 


0.2884 (46) 


1.305(21) 



Table B.8: Zo(go, L/a) at scale L = 1.436ro and Z^ GI (<?o) for both O12 and O44 and for = (1,0,0). Results are 
presented for the xSF with standard Wilson fermions and the SF with improved and standard Wilson fermions. 
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